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A major issue in the understanding of cuprate superconductors is the nature 
of the metallic state from which high temperature superconductivity emerges. 
Central to this issue is the pseudogap region of the doping-temperature phase di- 
agram that extends from room temperature to the superconducting transition. 
Although polarized neutron scattering studies hint at magnetic order associ- 
ated with the pseudogap, there is no clear thermodynamic evidence for a phase 
boundary. Such evidence has a straightforward physical interpretation, however, 
it is difficult to obtain over a temperature range wide enough to encompass both 
the pseudogap and superconducting phases. We address this by measuring the 
elastic response of detwinned single crystals, an underdoped YBCO6.60 with su- 
perconducting transition at T c = 61. 6K and a slightly overdoped YBCO6.98 with 
T c = 88. OK. We observe a discontinuity in the elastic moduli across the super- 
conducting transition. Its magnitude requires that pair formation is coincident 
with superconducting coherence (the onset of the Meissner effect). For both 
crystals the elastic response reveals a phase transition at the pseudogap bound- 
ary. In slightly overdoped YBCO that transition is 20K below T c , extending the 
pseudogap phase boundary inside the superconducting dome. This supports a 
description of the metallic state in cuprates where a pseudogap phase boundary 
evolves into a quantum critical point masked by the superconducting dome. 

Introduction. Resonant ultrasound spectroscopy (RUS) measures the frequencies f n of 
the mechanical resonances of a macroscopic specimen. 1,2 Using the redundant information 
contained in the frequency shifts A/ n (T), which in the measurements reported here are 
proportional to linear combinations of the small changes in elastic moduli, RUS can mea- 
sure the temperature evolution of elastic moduli Acy(T) with sensitivity exceeding 1 part 
per million. Like heat capacity, thermal expansion, and a very few other thermodynamic 
susceptibilities, elastic moduli are sensitive to all changes in the thermodynamic state of a 
system. If the pseudogap is a distinct phase like the superconducting transition, it must be 
accompanied by a thermodynamic signature evident in the elastic moduli. Unlike (scalar) 
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heat capacity, elastic moduli have several components, thus contain information about the 
symmetry of a phase transformation. The universal behavior observed in the metallic state 
at higher temperatures 3-7 does not extend below the pseudogap phase boundary, T*. Evi- 
dence for time reversal and inversion symmetry breaking 8-10 in the pseudogap phase suggests 
a continuous (second order) transition as developed in Ref. 11. 

A strength of ultrasonic measurements is that the thermodynamic information is accom- 
panied by transport information contained in the ultrasonic attenuation 12 , which unlike 
electrical conductivity is not hidden by superconductivity. For resonance methods, attenu- 
ation is measured by the inverse quality factor, proportional to the width of the resonance. 
Such measurements on millimeter-sized detwinned YBCO crystals face a range of experi- 
mental challenges. These include vibration isolation over a broad temperature range, stable 
specimen positioning without adhesives, use of minimal ultrasonic power while maintain- 
ing required signal-to-noise ratio, and temperature control under conditions where thermal 
contact to the crystal is weak (Figure 1; SI contains experimental details). 

Ultrasound and high temperature superconductivity in cuprates. Elastic moduli, 
second derivatives of the free energy, must change discontinuously (jump) across a second 
order phase transition, softening upon cooling into a symmetry-broken phase. The softening 
originates in the difference between the "fast" and "slow" elastic response in the symmetry- 
broken phase, —Acij = c^- — cfj. Analogous to adiabatic and isothermal measurements, here 
the term "fast" moduli, c?-, corresponds to the response at fixed order parameter rj, whereas 
"slow" moduli, cfj, correspond to moduli where the order parameter is at equilibrium, i.e. 
<p = dF/drj, the order parameter restoring force is zero. The order parameter has a finite, 
temperature dependent relaxation time r(T) which lengthens close to a transition (critical 
slowing down) 16 . This means that in a measurement at finite frequency, u, the softening 
extends over a finite temperature range in which the system crosses from l/r(T) <C u to 
1/t(T) ^> u. Causality requires that the dynamic origin of the broadening of the disconti- 
nuity in the frequency shift is accompanied by increased attenuation over the same range of 
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FIG. 1. (a) Mechanical resonance frequency shift versus temperature for underdoped YBCO6.60) 
T c = Ql.QK (hole doping p = 0.108±0.005 13 ), green, and slightly overdoped YBC0 6 . 9 8 (T e = 88K, 
hole doping p = 0.188), orange. Crystals are approximately 200/im thick, and one mm square. The 
two curves are offset vertically for clarity. In YBCO (as in most crystals) the largest component (up 
to 5%) of the temperature dependence of elastic moduli upon cooling from room temperature to 
zero is associated with the anharmonicity of the lattice 14 . The anharmonicity-driven frequency shift 
is smooth and monotonic 15 . (b) a typical resonance taken with the low acoustic power required, 
in-phase (red), quadrature (blue) components, (c) transducer-sample geometry, (d) balsa- wood 
probe used in the measurement. 



temperature (Kramers-Kronig in SI). 



Compressional elastic moduli (crystallographic scalars 17 ) enter linearly in the free energy 
and therefore can change discontinuously at a phase transition. Volume preserving moduli 
(shear) can be discontinuous only in their derivatives (break in slope). The magnitude of the 
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FIG. 2. (a,c) Resonance frequency with anharmonic background subtracted (a) and attenuation 
(c) across the superconducting transition in the YBCO6.60 crystal. Inset in (c) shows frequency 
shift across T c . 50mK intervals were used on cooling (blue), 70mK on heating (red) in a continuous 
temperature-controlled linear-in-time temperature sweep (SI for details). A small horizontal offset 
between heating and cooling, ~ 0.2K, clearly visible at the superconducting transition, is caused 
by imperfect thermal contact between thermometer and crystal. The origin of the vertical offset is 
unclear but not surprising for such high resolution scans. (b,d) Resonance frequency shift with an- 
harmonic background subtracted (b) and attenuation (d) in 90mK steps for the YBCOg.gs crystal. 
The maximum in attenuation in the YBCO"6.98 crystal is not associated with the superconducting 
transition as will be discussed in the text. Note that changes in elastic moduli clearly associated 
with the superconducting transition extend as far as 50K above and below T c . 



(dimensionless) step discontinuity, required at a the superconducting transition, is, within 
a factor of order unity, equal to (T C /T F ) 2 , where T F is a degeneracy temperature of the 
metallic state from which superconductivity emerges - for conventional metals the Fermi 
temperature. (SI contains details) 
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The temperature dependence of a given resonance /„ depends on changes in several elastic 
moduli Cij weighted with geometric factors afj, of order unity, which depend on crystal 
shape and the mode of vibration, Af n (T)/f n = a y(Acy(T)/cy). The lowest resonant 
frequencies are typically determined by shear elastic moduli alone, which in YBCO are 
about 1/3 that of compressional moduli 18 . The crystal shape determines the particular 
mix of moduli for each resonance. 1 The largest frequency shifts across the superconducting 
transition can be used to infer a lower bound on the thermodynamic effects at the phase 
transition. In underdoped YBCO6.60 we measured a shift Af / f w 10~ 4 , and for slightly 
overdoped YBCO, Af / f is an order of magnitude larger (Figure 2). Because there is 
no significant attenuation increase at the superconducting transition, the transition width, 
AT C = 0.5K, is not associated with dynamics 19 ' 20 and therefore must be caused by small 
variations in hole doping over the volume of the crystal. With Tp ~ 5000K and T c /Tp ~ 
1CT 2 , the magnitude of the observed elastic discontinuity at both dopings is consistent 
with the entire thermodynamic signature of superconductivity occurring within 0.5K of the 
transition. 

The phase transition at the pseudogap phase boundary. On cooling, the resonance 
frequency shift in YBCOq,qo reveals a break in slope at T* = 245K, a standard marker for a 
phase transition, Figure 3. This temperature corresponds well with the onset temperature of 
magnetic order reported in YBCO specimens of similar composition. 9 For phase transitions of 
the class indicated by neutron diffraction 9 ' 10 there is no requirement for a step discontinuity. 21 
The pseudogap temperature, T*, varies much more rapidly with oxygen composition than 
the superconducting transition temperature - near oxygen doping x = 6.60, dT*(p)/dp is 
about 10 times larger than dT c (p)/dp, Figure 3(d). This is consistent with the observed 
AT* = 5K width of the pseudogap transition, about 10 times that of the superconducting 
transition, (dT* / dT c )AT c . 

The phase transition at T* = 2A5K is accompanied by a strong increase in attenuation above 
T*. This ultrasonic attenuation 12 is caused by the slowing down of fluctuations, reaching a 
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maximum when the measurement frequency matches the characteristic timescale, ujt(T) = 1 
(details in SI). Frequently 16 attenuation reaches a maximum on the symmetry-broken side 
of a transition and spans a temperature range equal to the transition width. However, 
for YBCO (Figure 3) the attenuation maximum occurs at a temperature higher than the 
pseudogap transition, and spans a temperature range wider than the width of the transition. 
On heating, the pseudogap transition is broader than on cooling. The transition exhibits 
approximately 10K hysteresis (details in SI). The small discontinuities in frequency observed 
upon cooling in the pseudogap phase (Figure 3(a)) do not occur in similar sized well behaved 
metals in the same apparatus. Such variations are reminiscent of the response of domains, 
especially in that they depend on cooling rate and history and are absent on heating. It is 
possible that although the crystals are detwinned, a very small number of crystallographic 
twin boundaries could be present. 

For the slightly overdoped crystal, YBC06.98, the break in slope occurs at 67K. Extrapolating 
T*{p) from measurements in underdoped cuprates (Figure 3(d)), we conclude that this phase 
transition is the pseudogap phase boundary, extending T* 20K below the superconducting 
transition near optimal doping. This extends the doping range over which magnetic order 
exists to the vicinity of optimal doping. The width of the attenuation maximum associated 
with the pseudogap appears to scale with T*. The resonances exhibit a strong (up to 
hundred-fold) increase in ultrasonic attenuation over a narrower temperature range (60K and 
80K) than in the underdoped crystal. This reveals a correlation between the measurement 
frequency and the temperature of the attenuation maximum (Figure 4). Interpreting this 
frequency dependence as originating from ujt(T) = 1 reveals the divergence of r(T) as the 
pseudogap phase transition is approached. In Figure 4(e) the extrapolated value of l/r(T) 
vanishes close to 65K. Like the underdoped specimen, the strong increase in ultrasonic 
attenuation is observed above the pseudogap phase boundary. Note that the frequency shift 
is upward on cooling to the symmetry broken phase (SI for more discussion). 

The resonance measurements presented here provide quantitative bounds on the dis- 
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FIG. 3. Resonance frequency shift (a) and linewidth (b) vs temperature in YBCO6.60 f° r several 
modes. For clarity the frequency shift curves are offset vertically. A clear increase in ultrasonic 
attenuation is observed close to the onset temperature of Kerr rotation 22 . The many different 
temperature dependences in (a) can be decomposed into three linearly independent components 
(d) (SI has details). It is required that there be no more than 6 such components (tetragonal 
symmetry) and, likely, fewer because some of the six independent elastic moduli will have similar 
temperature dependences. At least one should not show anything at T* , as is observed, (c) Phase 
diagram of hole-doped cuprates. Neutron diffraction 9 (blue squares), magnetic susceptibility 23 
(triangles), Kerr rotation 22 (circles), Resonant Ultrasound (stars). The error in determining T* 
from neutron diffraction is about 50K. The error in determining T* from RUS measurements is the 
apparent width of the transition, not the intrinsic precision of RUS (frequency errors are negligible, 
temperature errors are less than 20mK). 

continuity in the elastic moduli across the superconducting transitions, and rule against the 
thermodynamics of pairing occurring at temperatures above T c in both underdoped and 
slightly overdoped YBCO. These observations also support a description of the metallic 
state and high temperature superconductivity in cuprates where a pseudogap phase bound- 
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FIG. 4. RUS measurements across the pseudogap phase boundary in slightly overdoped YBCO. 
Ultrasonic attenuation (a) and resonance shift vs temperature (b) for three resonance peaks showing 
strong (up to hundred-fold) increase in ultrasonic attenuation. There is a weak maximum near 
120K that is frequency dependent, (b) curves offset vertically, (c) logarithm of resonance width 
for several modes, with vertical offset. Gray: evolution of the attenuation maximum with mode 
frequency, (d) independent linear components of the temperature dependence of several resonance 
peaks. The green curve shows a break in slope at T*=67K, the red curve shows only a smooth 
background also present in green and purple, and the purple curve shows the superconducting 
discontinuity, (e) from the plot of attenuation versus temperature (c), we extract the temperature 
of the attenuation maximum, and the frequency at which this occurs (noting that ujt(T) = 1 at 
that frequency) and plot l/r(T) versus that temperature (e). 

ary evolves into a quantum critical point masked by superconductivity. 
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